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Abstract 

Supersonic expansion is a very powerful tool to produce an atomic beam with a well defined 

,^ ■ velocity and, by seeding a test gas in such an expansion, the energy of the test gas can be transferred, 

'p^. at least partially, to the very-low-temperature carrier gas. The case usually studied is the one of a 

heavy gas seeded in a light carrier gas and, in this case, the parallel temperature of the seeded gas 

is always larger than the one of the carrier gas. In the present paper, we study the opposite case 

d ' which has received less attention: when a light gas is seeded in a heavier carrier gas, the parallel 

c/3 ■ temperature can be substantially lower for the seeded gas than for the carrier gas. This effect has 

• i— I . been first observed by Campargue and coworkers in 2000, in the case of atomic oxygen seeded in 

J>-^, argon. In the present paper, we develop a theoretical analysis of this effect, in the high dilution 

i-^ \ limit, and we compare our theoretical results to several experimental observations, including a set 

of measurements we have made on a beam of lithium seeded in argon. The agreement between 
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theory and experiments is good. 



I. INTRODUCTION 

We have recently observed that, in a supersonic beam of argon seeded with hthium, the 
parallel temperature of lithium is roughly one third of the parallel temperature of argon [l| . 
A similar effect was reported in 2000 by Campargue and co-workers J2] on a beam of atomic 
oxygen seeded in argon: the oxygen parallel temperature was close to half the argon parallel 
temperature. This result was found to be in agreement with a calculation developed by the 
same authors. 

This effect is surprising because one would naively expect that the parallel temperature 
of the seeded gas cannot be lower than the parallel temperature of the carrier gas: the 
expansion of the carrier gas is the source of the cooling effect and the energy of the seeded 
gas is transferred by collisions to the carrier gas, but the transfer should not be complete. 
This view is too naive as shown below, but it is supported by all the experiments involving a 
heavy gas is seeded in a light carrier gas: the parallel temperature of the seeded gas always 
exceeds the one of the carrier gas and the ratio of these two temperatures increases steadily 
with the mass ratio, from 2.2 — 2.5 for an He-Ar mixture, up to 4 — 4.5 for He-Xe mixture, 
and reaching 13 — 15 for H2-Xe mixture, following results due to Campargue and co-workers 
J3|, y|. At the same time, the mean velocities of the two components are different and this 
difference, the velocity shp effect p, |6| , increases also with the mass ratio. 

The theory giving the terminal temperatures in supersonic beams of a pure monoatomic 
gas was developed by many authors including Anderson and Fenn J7| , Hamel and Willis J8[ , 
Knuth and Fisher [9|, Miller and Andres [6|, Toennies and Winkelmann [lO| and by Beijerinck 
and Verster [11]. The extension to the case of binary mixtures was made by Cooper and 
Bienkowski |lj| , by Anderson and coworkers la |lj] , by Chesneau and Campargue [17[ (the 
subject has been reviewed by D.R. Miller (l3|). Because a large majority of experiments 
corresponds to a heavy species seeded in a light carrier gas, the calculations were done in 
this case (reference J2| being an exception). These calculations established that the terminal 
parallel temperatures of the two species are not equal jlj| and were able to explain the 
observed ratio of parallel temperatures of the seeded and carrier gas [15|, [ly, |l7||- In the 



present paper, we apply this theory to the case of mixtures of monoatomic gases. We need 
several approximations to make an analytic theory and the most important one consists in 
neglecting the velocity slip effect. This approximation is good when the seeded gas atomic 
mass 7712 is smaller than or comparable to the one of the carrier gas mi. We are able to 
calculate the ratio of parallel temperatures of the seeded and carrier gases, as a function of 
their interaction potentials and atomic masses. 

This theory provides a physical explanation of the ultra cooling effect: during a supersonic 
expansion of a pure gas, the basic phenomenon is the geometrical cooling effect of the 
perpendicular temperature, as illustrated by figure 1 of reference jlO|, and this cooling effect 
is transferred by collisions to the parallel temperature. When a gas is seeded in such an 
expansion, with a high dilution, the parallel temperature of the seeded gas is coupled to its 
own perpendicular temperature as well as to the parallel and perpendicular temperatures 
of the carrier gas. At the end of the expansion, the perpendicular temperatures of the 
seeded and carrier gases keep on decreasing and become considerably lower than the parallel 
temperature of the carrier gas which is frozen. Therefore, the terminal parallel temperature 
of the seeded gas can become lower than the same quantity for the carrier gas, provided that 
the collisional coupling of the temperatures acts for a longer time for the seeded gas than 
for the carrier gas. We must compare the efficiency of collisional exchange of kinetic energy 



during the various possible collisions. Two effects govern this comparison: i) the relative 
ranges of the seeded gas-carrier gas interaction and the interaction between two atoms of 
the carrier gas; ii) the mass ratio as it plays a very important role in the way energy is 
redistributed during a collision. Among these two effects, the mass ratio usually dominates 
and we predict that the parallel temperature of the seeded gas can be considerably lower 
than the parallel temperature of the carrier gas when 1712 <^ rrii. 

The present paper is organized as follows. We first recall the theory of the terminal 
temperatures in a pure gas expansion, thus introducing the needed ideas, notations and 
equations. Then, we generalize this calculation to the case of a gas mixture, in the limit of a 
high dilution. In the final part, we compare the available experiments with our calculations. 

II. TERMINAL PARALLEL TEMPERATURE IN SUPERSONIC EXPANSION 
OF A PURE MONOATOMIC GAS 

In this part, we briefly recall the theory following the 1977 paper of Toennies and Winkel- 
mann llOll fnoted below TW) and we introduce the simplifications used by Beijerinck and 
Versternni in 1981 (noted below BV). 

A. Equations describing the cooling effect during a supersonic expansion 

The starting point is the Boltzmann equation in the steady state regime: 

vgradKr)/(r,v)]=[^M'j (1) 

V / coll 

n{r) is the gas density and /(r, v) the normalized velocity distribution. Following TW, the 
supersonic expansion is approximately described, near its axis, by a spherically symmetric 
flow. The velocity distribution is assumed to remain Maxwellian, with different parallel and 
perpendicular temperatures T\\ and T±: 
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where m is the atomic mass and u is the local hydrodynamic velocity. Then, TW express the 
Boltzmann equation in spherical coordinates and apply the method of moments to obtain 
four differential equations coupling the density n, the hydrodynamic velocity u and the 
temperatures Ty and T± (see equations collected in table I of TW). After some algebra and 
neglecting small terms proportional to kBT\\/{mu^) = 1/(25*1?) where S\\ is the parallel speed 



ratio S\\ = u/ J2kT\\/m (a good approximation as 6*11 becomes large in most supersonic 
expansions), one obtains two differential equations coupling the parallel and perpendicular 
temperatures: 
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where, following BV notation, we use z to measure the distance from the nozzle. JF is a 
collision term given by: 

^ = ^ J 9^AEfiy^,)fiv,)d\,d\,dn (5) 

In equation (0), Vi and V2 are the atom velocities before the collision and g = vi — V2 is 
their relative velocity of modulus g, while da{g)/dQ is the differential cross-section. AE 
is the energy transferred during one collision from the parallel degree of freedom to the 
perpendicular ones for the two atoms. After averaging over the azimuth describing the 
direction of the final relative velocity around the initial relative velocity, AE is given by : 

{AE) = j[gl-2gf^ [l-cos^x] (6) 

where x is the deflection angle. 



B. Simplification of these equations 

BV have shown that the coupling term JF is well approximated by a linear function of 
(2]| - T^) given by J^ ^ A{z) (t,, - T^) /2 with A{z) = 16n{z)n'-^''^\Tm) / {l5uoc) ■ Here, 

f2^''*-'(T) is a thermal average of the collision cross-section Q^''\ both defined in reference |lq| 
(see also Appendix A). T^ = (Ty + 2Tj_)/3 is the weighted mean of the parallel and per- 
pendicular temperatures. The hydrodynamic velocity u, in equation (0), has been replaced 



by its terminal value Mqo = J^kBTQ/m in the large S\\oo limit. Neglecting quantum effects, 
assuming a 12-6 Lennard- Jones potential, the il'^^'^)(T) integral is given by: 

n<-.,T,.2.99(^)"(£^)" (7) 

where Cq is the coefficient of the attractive term of the potential. This formula, established 
by BV, is valid when fc^T is small with respect to the potential well depth e (see appendix 
A). n{z) is related to the source density uq and temperature Tq by n{z) ^ I /{uooZ^) with 
the intensity I of the supersonic beam given by / = UQUooZ^^r with z^ef = 0.403 x d {d is the 
nozzle diameter): 
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with A{z) given by: 
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A,„. 3.189 x-f-- (^) (10) 



C. Scaling and integration of these equations 

When z is small, the density ni^z) and the collision term JF are both large and the two 
temperatures remain equal, Ty = T^ = T^ which verifies: 

~dz ~ 37 ^ ' 

giving Tm oc z"'^!'^ . BV simplified the equations ()8I9|) by introducing reduced temperatures, 
obtained by dividing the temperatures Ty, T^, T^, by Tq and a reduced distance Zr = zjzref- 

5^ = -S% (T||. - T^.) (12) 

dTxr 2Ti^ T^/.^ / N 

All the source parameters are condensed in the quantity S = 0.813 xrioc? {C^lksT^ , which 
can be eliminated by a further scaling due to BV, Zr = ^S^/^^ and Txr = rxS"^^/^^ (with 
X = II or ±), thus providing universal equations: 

dT± _ 2t^ rV^(r|| - r^) . . 

^ " C^ ^ ^ 

These equations have been integrated numerically by BV and, using MATLAB, we have 
reproduced their calculation (see figure 1 below). When ( is large, t± oc (^^ while ry tends 
toward a limit, ryoo ~ 1.15, from which we can express the final parallel temperature as a 
function of source parameters: 

T||oo/To = 1.151 X S-12/11 (Ig) 

D. Results and tests in the case of argon expansion 

The tradition is to give the terminal value 5*1100 of the parallel speed ratio rather than 
the terminal parallel temperature, with the relation Tyoo ~ 2.5Tq/S'u^. We have written the 
results of TW and BV in the same form to facilitate their comparison: 
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TW obtained A = 1.413 and 6 = 0.53 by fitting the results of numerical integration of their 
equations written without approximations. Following the procedure we have just recalled, 
BV obtained A = 1.313 and 6 = 0.545. In the range of practical interest, when 5 < S\\oo < 50, 
these two theoretical formula never differ by more than 4% and agree when S\\oo ~ 18.9. 

From a fit of experimental S\\oo values for argon, BV obtained semi-empirical values of 
the parameters A = 1.782 and 6 = 0.495, using CeiAr - Ar)/kB = 4.45 x 10"^^ K.m*^. We 
have verified that this semi-empirical formula represents very well the ^Hoo values for argon 
measured by H. D. Meyer [191 and we estimate the error bar on 5*1100 near ±7%. 



III. GENERALIZATION TO THE CASE OF A MIXTURE OF TWO 
MONOATOMIC GASES 

A. Approximations used to solve the Boltzmann equation 

We have two Boltzmann equations, one per species, noted i = 1 for the carrier gas and 
i = 2 for the seeded gas and the colhsion terms have two parts corresponding to the two 
colhsion pairs: 

V, ■ grad K(r)/,(r, v,)] = S,=i,2 (^^) ihj) (18) 

V / coll 

The densities /j(r, Vj) are expressed by equation (0), with the atomic masses mj, the tem- 
peratures T||j and T^j. We assume that the hydrodynamic velocity u is the same for both 
species. When a heavy species is seeded in a light gas, the velocity slip effect J5|, |6| is not 
negligible and this approximation would be bad, but we are interested in the opposite case, 
when the atomic mass of the seeded gas is smaller than or comparable to the atomic mass 
of the carrier gas. When the expansion is well in the supersonic regime, the difference of the 
mean velocities is negligible with respect to the thermal velocity in the moving frame (see 
below). 

We consider the high dilution limit, when the seeded gas density ^2 is considerably smaller 
than the carrier gas density rii. Because 712 <^ rii, in the Boltzmann equation for the carrier 
gas distribution function /i, we can neglect the 1 — 2 collision term. The expansion of the 
carrier gas is not modified by the seeded gas and the equations written above can be used 
to calculate the parallel and perpendicular temperatures of the carrier gas. We must simply 
introduce in A{z) the relevant collision integral fl\{ (T), where the indices designate the 
colliding atom pair. In the Boltzmann equation for the seeded gas /2, we consider only the 
effect of 1 — 2 collisions, i.e. collisions with atoms of the carrier gas and we neglect the 
collisions involving two atoms of species 2. Then, we get the following equations for the 
temperatures of species 2: 

2J^II,2 (19) 

+ -^±2 (20) 
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with the coUisional energy transfer terms given by: 

ni f ^dai^2{g) 

where Ai?x,2 {X = || or ±) measures the parallel or perpendicular energy gained by atom 2 
during the collision with atom 1. We can express A£^x,2 with the parallel and perpendicular 
components of the center of mass velocity and of the relative velocity. However, when the four 
temperatures are not equal, the product /i(vi)/2(v2) has not a simple form when expressed 
with these velocities and we have not been able to calculate exactly the resulting integrals. 
Therefore, to get an analytic result for J^x,2i we have used the following approximation. 
When the collision energy is small, the cross-sections Qi2{g) behave like g'"^^^ and the 



products gQi^2i9) vs-^y slowly with g, like g^^^. We take these products gdai^2{g)/d^ out of 
the integrals over vi and V2 to get: 
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The calculation of this integral is described in appendix B. 
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B. Coupled equations describing the temperatures of the two gases 



The coupled equations for the carrier gas i = 1 have been established in part II (equa- 
tions (jHI) and 0) where Tx must be replaced by Txi- In the differential equations for the 
temperatures of the seeded gas « = 2, we have quantities similar to K{z) but involving the 
Vti2 [Tm) integrals with / = 1 and 2. By introducing two dimensionless ratios ps and po, 



(2,2), 



)(2,2) 



we can express these quantities as a function of Vl\^{ (Tm)- Ps is the ratio of ^Ij collision 
integrals differing by the species of the second collision partner and its value is deduced from 
equation ((7j): 
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ps depends slowly on the Cq ratio and more rapidly on the mass ratio, po is the ratio of 
angle-averaged cross-sections of orders / = 1 and / = 2, defined by po = ^[^2 /^i,2 = 1-32 
(see Appendix A and our previous paper [l|). 
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A{z) being given by equation (fTnjl and Tx,av being defined by: 



(25) 



Tx,av = PTxi + aTx2 (26) 

with a = TTLi/M, (3 = m2/M, M = mi + 1712 and p = mii7i2/{mi + 1712). One should remark 
that in Tx,av, the weight of Txi is /3 = 1712 /M and the weight of Tx2 is a = mi/M. We 
proceed as in part I.C and we get two equations verified by rp and t±2'- 
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A simple test of the coherence of our calculations is to consider that the two species have 
the same masses {mi = 7712) and the same collision cross-sections {ps = 1). We then expect 
that the parallel and perpendicular temperatures are independent of the species (rxi = tx2 
for all ( values) and this property is well verified. 

Using MATLAB, we have integrated the equations ()27p28j) and the results are represented 
in figure 1. At the end of the expansion, the parallel and perpendicular temperatures 
of the seeded gas are intermediate between the same quantities for the carrier gas. In 
particular, the terminal value of the parallel temperature ratio T||2, 00/^111, 00 = ^|2,oo/^|i,oo 
can be substantially lower than 1. The absolute minimum value for this temperature ratio is 
reached when the right-hand side of equation (j^3j) vanishes. Assuming that the perpendicular 
temperatures are both negligible, we get the minimum possible value of the ratio T'||2,oo/7]|i,oo: 

Min C^^] = ^^^° (29) 

V^|i,oo/ mi + 2m2Po 

This limiting value would be reached if the ratio ps tends toward infinity. We have plotted 
in figure 2 the terminal value of the ratio T||2,oo/r|ii,oo for various values of the ratio ps, as a 
function of the mass ratio 7712/777,1. When the mass ratio 1712/ mi is very small, the parallel 
temperature of the seeded gas can be considerably smaller than the one of the carrier gas. 
We can make some surprising predictions: 

• In the case of molecular hydrogen H2 seeded in argon, using to calculate the Ar-H2 Cq 
coefficient the combination rule and the data of Kramer and Herschbach |2fl] , we get 
Ps = 2.48 from which we predict a temperature ratio 2^112, oo/7'||i,cx) = 0.16, a quite large 
effect! 

• for a beam of potassium seeded in argon, the two atomic masses are almost equal and, 
with the Cq from reference |26|, we get ps = 1.67. We predict a temperature ratio 
^||2,oo/^|i,oo = 0.77, already smaller than 1. 

• for a beam of sodium seeded in neon, we calculate in the same way ps = 1.85, (the 
large Cq ratio overcompensates the effect of the mass ratio 7712/ m,i = 1.14) and we 
predict a temperature ratio T'||2,oo/7'||i,oo = 0.77, also smaller than 1. This is not a 
large effect, but this result, with T||2,oo < ^||i,oo while 7712 > tui, is contrary to the 
common belief recalled in the introduction. 



C. Discussion 

In addition to the approximations already done by Beijerinck and Verster, we have ne- 
glected the velocity slip, which has been studied in particular by Anderson [5| and by Miller 
and Andres pj. Equations (2.37,2.38) of Miller [13|) give the velocity difference M2,oo — ^i,cx), 
which can be written with our notations and in the high dilution limit: 

8 
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physically important quantity is the ratio of this velocity difference (u2,oo — mi^oo) divided 



The 0,1 2 integral, which appears in these equations, is evaluated in Appendix A. The 



by the parallel thermal velocity of atom 1, namely J2kBT\\i/mi. Using equation ()17|) in the 
TW form and approximating the 1.07 exponent by 2 x 0.53, we get: 



M2,oo - Ml,oo 0.84 
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(31) 



If the terminal parallel speed ratio S'nioo is large, the velocity slip is a small fraction of the 
parallel thermal velocity of atom 1, provided that the masses irti and m2 are not extremely 
different. In this case, neglecting the velocity slip is an excellent approximation. 

In equation (PT|) . we have also treated the products gQi^2i9) &s constant. Takahashi and 
Teshima were able to calculate numerically these integrals when they studied the case of a 
heavy gas seeded in a light gas [16|. Using their technique, it would be possible to test the 
accuracy of our approximation. 

Monte Carlo simulation can also be used to simulate flows of gases and gas mixtures. The 
Direct Simulation Monte Carlo method described by G. A. Bird in his book ■22.1 is the best 



method and it has been used very early to simulate supersonic expansions |l5l . l23l |. This 



method requires much computation in the high dilution case J21| and P. A. Skovorodko has 
used the Test Particle Monte Carlo method to simulate the expansion of gas mixtures. This 
method also involves an approximation, the use of Maxwell type interaction potential. He 
has simulated our lithium seeded in argon expansion and his results support our calculations, 
at least at the qualitative level |2J|. In his simulations, P. A. Skovorodko has found that the 
distribution of the perpendicular velocity of lithium differs very much from a Maxwellian 
distribution and is almost perfectly exponential. The fact that the velocity distribution of the 
perpendicular degree of freedom cannot be Maxwellian was analyzed in detail by Beijerinck 
and co-workers who suggested to call blistering this effect [ll|, |25| . This effect cannot be 



taken into account by the theory described in the present paper as the assumption of an 
elliptic Maxwellian velocity distribution is the starting point of the solution of Boltzmann 
equation. 

IV. COMPARISON WITH EXPERIMENTAL RESULTS 

In this part, we first describe briefly the previous experimental cases in which the light 
gas parallel temperature has been found lower than the carrier gas parallel temperature. 
Then, we describe our experiment with a beam of lithium seeded in argon and a new set of 
temperature measurements. 

A. Sodium in seeded in argon 

A beam of sodium seeded in argon was built by D. Pritchard and coworkers ^33] with a 
source temperature near 1000 K, an argon pressure up to po = 3 bars and a nozzle diameter 



d = 70 /iin. In a particular experiment j3^, the mean velocity was measured, u = 1040 ± 2 
m/s corresponding to Tq = 1039 K, and the rms velocity width, deduced from the atomic 
diffraction pattern, was found equal to 3.7 ± 0.4% , corresponding to T||2oo = 4.1 ± 0.9 
K. Assuming that this experiment was made with the largest pressure quoted in 33[, we 



calculate the argon parallel temperature Tyi = 7.7 ± 1.1 K, from which we deduce the 
parallel temperature ratio T||2oo/^|ioo = 0.53 ± 0.19. This value is in good agreement with 
our theoretical result is T||2oo/7'||ioo = 0.66, deduced from the value of ps = 1.67 (obtained 
with the Cq values from reference ;2Q]). 

B. Atomic oxygen seeded in argon 

In 2000, Campargue and co-workers J2| operated and characterized a beam of atomic 
oxygen seeded in argon produced by a laser plasma source. The oxygen parallel temperature 
was roughly half the argon parallel temperature (see figure 10 of reference 2]). The numerical 
model described in this paper supports the experimental results. Let us compare this result 
with our theory. Usin g th e combination rule of Kramer and Herschbach J20| and data for 
oxygen from reference [23], we calculate the oxygen-argon Cq coefficient Cq{0 — Ar) = 34 
atomic units and ps = 1.06. Our theory then predicts a temperature ratio ^112,00/^111,00 = 
0.76, substantially larger than the experimental value near 0.5. This discrepancy is probably 
due to the fact that, because of the very high source temperature close to 10"^ K, the 
terminal parallel temperatures are considerably higher in this case than for usual beam 
experiments, in the 70 — 130 K range for oxygen and 220 — 150 K range for argon. For 
such terminal temperatures, we cannot expect the low temperature approximations of the 
Q(^'^)(T) integrals to be valid. 

C. Lithium seeded in argon 

1. The beam source 



Our beam j28|] is inspired by the design used by Broyer, Dugourd and co-workers to 
produce lithium clusters [29j. The temperatures used in our experiment are usually equal 
to 973 K for the back part of the oven, corresponding to a lithium vapor pressure of 0.55 
millibar and to 1073 K for the front part. The nozzle is a hole of 200 pm diameter drilled 
in a stainless steel 0.3 mm thick wall. The argon gas (from Air Liquide, 99.999 % stated 
purity) is further purified by a purifying cartridge also from Air Liquide and its pressure 
can be varied from 150 to 800 millibar, limited by the throughput of our oil diffusion pump 
(Varian VHS400 with a 8000 1/s pumping speed) backed by a Leybold D65B roughing pump 
(65 m^/hour pumping speed). For a source pressure po = 300 millibar, the pressure at the 
diffusion pump is 8 x 10~^ millibar. We use a skimmer from Beam Dynamics with 0.97 mm 
aperture at a 20 mm distance from the nozzle. After the skimmer, the lithium beam is in 
a separate vacuum tank pumped by an oil diffusion pump (Varian VHS6 with a 2400 1/s 
pumping speed) fitted with a water cooled baffle. Under beam operation, when po = 300 
millibar, the pressure in this chamber is 3 x 10~^ millibar. 
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2. Doppler measurement of the parallel and perpendicular velocity distribution 

In the center of the second vacuum tank, i.e. 225 mm after the skimmer, the hthium 
beam is crossed by two laser beams A and B. The angle between the atomic beam and the 
laser beams are 6a = 47.9 ± 0.5° and 6b ~ 90°. The first order Doppler effect is sensitive 
only to the projection Vp of the velocity on the laser beam axis. For a laser beam making the 
angle 6 with the axis of the atomic beam, the distribution of Vp is deduced from equation 

/ \ 1/2 

^("-) = ( 2^^£b) ) ^"P 



m[Vp — u cos ( 



2ksn6) J ^^^^ 

with T{6) = T|| cos^ 6 + T± sin^ 6. The distribution f{vp) is centered at Vp = u cos 6 and its 
width is characterized by a weighted mean of the parallel and perpendicular temperatures. 
The fluorescence intensity as a function of the laser frequency reflects the velocity distribu- 
tion if the natural width of the excited transition is negligible with respect to the Doppler 
width and if saturation broadening of the transition as well as laser frequency jitter are both 
neghgible. 

The laser excites successively the hyperfine components of the ^5*1/2 - ^-P3/2 resonance 
transition of lithium at 671 nm J30|. The very small hyperfine splittings of the upper state 
can be neglected. The ground state has two hyperfine components F = 1 and F = 2, 
with a splitting equal to 803.5 MHz. The natural width of the transition is r/27r = 5.87 
MHz J3l| and we use laser power density of the order of 10~^ mW/cm^, corresponding to a 
saturation parameter s ~ 5 x 10~^ and a negligible broadening of the excitation line. Our 
single frequency cw dye laser pumped by an argon ion laser has a linewidth of the order of 
1 MHz, thanks to the Hansch-Couillaud [3^ frequency stabilization technique. 

Figure 3 shows such a laser induced fluorescence spectrum: for beam A, the Doppler full 
widths are of the order of 200 MHz and the natural width as well as the laser linewidth are 
negligible. During our experiment, we have also recorded a saturated absorption spectrum 
in a heat pipe oven: this signal provides Doppler-free peaks which are useful if one wants to 
measure accurately the mean velocity u of the atomic beam, because, in some experiments, 
the laser beam B is not exactly perpendicular to the atomic beam. 

3. Experimental results 

From the analysis of the fluorescence signals, we can deduce the beam mean velocity 
u = 1010 ± 10 m/s, slightly less than the value deduced from the source temperature. 



■^oo = y5kBTo/m = 1056 m/s. The width of the B peaks is related to the perpendicular 
temperature of the beam but we think that the perpendicular temperature T± thus deduced 
is overestimated because our detector is not observing only one streamline. The width 
of the A peaks provide the weighted mean T{6) of the temperatures, from which we can 
deduce the parallel temperature. As we think that our measurement may overestimate 
the perpendicular temperature, we have neglected the contribution of the perpendicular 
temperature when extracting the parallel temperature from T{6), i.e. we have used: 

T\\=T{6)/cos^6 (33) 

so that we get an upper limit of the parallel temperature Ty. In table I, we have also given 
the fluorescence signal intensity which gives an idea of the intensity of the lithium beam. 

11 



Source pressure (mbar) 


TpiK) 


T±2imK) 


Fluorescence signal(a.u.) 


200 


10.1 


493 


8.5 


267 


7.6 
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7.9 


333 
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7.0 


400 


6.0 
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467 


6.1 


546 


6.6 


534 


6.0 


593 


5.0 


600 


5.3 


550 


5.6 



This intensity decreases when the argon pressure increases. This behavior is probably due 
to the interaction of the beam with the residual gas near the skimmer. 

The measured parallel temperature is plotted as a function of argon pressure in figure 4. 
The full curve is a plot of our results with the ratio Tp/Tyi taken equal to its predicted value 
0.38. The dashed curve is obtained by fitting the value of this ratio Tp/Tyi to the four data 
points with an argon pressures smaller than 400 millibar: the fitted value is equal to 0.31. 
T||i is related to the source parameters by equation (fT7|). using the semi-empirical coefficients 
of BV. For larger pressures, the parallel temperature remains constant while theory predicts 
that it should still decrease. We think that this discrepancy may be due to the interaction 
of the beam with the residual gas discussed above. 

Another explanation could be cluster formation in the expansion. Following the results of 
Hagena ^] further discussed by Beijerinck and Verster ll|, for an argon beam, a safe upper 
limit to prevent cluster formation is to keep the quantity pQcf'^^Tf^^'"^ < 4.9 x 10^^ mbar ni°'^^ 
K~^'^. At our largest source pressure po = 600 mbar, the quantity podP'^^Tf^'^'^ = 3.6 x 10~^ 
mbar rn^*^^ K~^'^ is well below this upper limit. As the LiAr molecule is considerably less 
bound than the Ar2 molecule |.39- .40.] and as the lithium concentration is very small, we 
think that the threshold for the formation of mixed clusters is also not reached. 



V. SUMMARY 



The knowledge concerning the terminal parallel temperatures produced by supersonic 
expansions of a mixture of monoatomic gases has been developed for many years and the 
case of heavy atom seeded in a light carrier gas, which has been subject of almost all the 
studies, is extremely well understood. The opposite case (a light gas seeded in an heavier 
carrier gas) had been rarely studied, an exception being the work of Campargue and co- 
workers [2| in 2000. This opposite case might seem of limited interest but some experiments 
(for instance, atom interferometry) require a rather slow atomic beam with a small velocity 
dispersion: in this case, seeding a light gas in an heavier carrier gas is the most natural 
solution. 

In the present paper, we have developed a theoretical analysis by making a slight extension 
of the theoretical works of Toennies and Winkelmann jlO] and of Beijerinck and Verster jll| . 
With some new approximations, we obtain an analytic treatment of the problem, valid in the 
high dilution limit. We thus rationalize the surprising experimental result we had observed 
after Campargue and co-workers [2]: the parallel temperature of the seeded gas can be lower 
than the same quantity for the carrier gas. We have found that this effect depends on two 
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quantities, the ratio ps of cross-sections for the two gases, defined by equation (|23|). and the 
ratio of the atomic masses 7712/1711. Finally, the physical basis of this ultra cooling effect is 
not mysterious and we have explained it in the introduction. 

In a final part, we have discussed the available experimental evidences of this effect, due 
to other groups J2, |33|, |3J| and we have described our experiment with a beam of lithium 
seeded in argon. We have presented a set of experimental data considerably larger than in 
our first study jl|. The agreement between our theoretical prediction and experiments can 
be considered as good in the cases of sodium and lithium seeded in argon. The agreement 
is less good for the argon-oxygen beam but we understand why it is so. 

We think that the ultra cooling effect described here is now well established and some 
further studies would be very interesting. First of all, it should be easy to observe this 
effect in several other experimental cases, some of them being discussed in part III.B: a 
particularly striking case would be to seed molecular hydrogen in argon or even heavier 
rare gases. From the theoretical point of view, it is interesting to investigate this effect by 
other methods and Monte Carlo simulation, which has been already applied to this problem 
by P. A. Skovorodko 2J|, is particularly well suited as the approximations needed to solve 



Boltzmann equation can be removed. 
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VII. APPENDIX A: WEIGHTED CROSS-SECTIONS; THEIR DEFINITION 
AND CALCULATION OF THEIR TEMPERATURE DEPENDENCE 

A. Definitions 

We consider the collision between two atoms of species i and j with a relative velocity 
g and a reduced mass fi = mimj/{7ni + 7nj). The definitions and many results are taken 
from the book Molecular theory of gases and liquids, by Hirschfelder, Curtiss and Bird [18i] . 
Starting from the differential cross-section daij{g)/dfl, the angle- weighted cross-sections 
Qf]{g) are defined by (equation 8.2-2 of |l3) : 



fl(p)=f^-^^(l-cos^ 



Ql;j(5)=y ^^(l-cos'x)rffi (34) 

where x is the deflection angle. In classical mechanics, the deflection angle x is a function 
of the impact parameter b and the relative velocity g. Qij{g) is then given by: 

Q%) = J {1- cos' x)2nbdb (35) 
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The thermal averages Q\f (T) are given by equation (8.2-3) of reference [18[: 



i^fiT) = ^kBT/{2ni^) /°°Qg(^)y^+3exp(-f )d7 (36) 



with 7^ = fig'^/{2kBT). Three circumstances contribute to make classical mechanics a good 
approximation in the experimental cases considered here: 

• the quantum character of the interaction between two identical atoms of mass m 
is measured by the quantum parameter i] = h? /{ma'^e) introduced by Stwalley and 
Nosanow J36[ (see also De Boer and Lunbeck [37[) where e is the potential well depth 
and a the core radius. A large r] value indicates a highly quantum behavior while 
a small one corresponds to a quasi-classical behavior. "^He dimer is highly quantum 
with r] ^ 0.18 ^36], while Ar dimer l40] with r^ pa 0.75 x 10"^ and '^Li - Ar with 
7] K, 6.6 X 10^'^ have a weak quantum character. 

• the differential cross-section presents several types of quantum effects. These effects 
are largely reduced by the angle average following equation (jM|) and further reduced 
by the thermal average following equation ()36|) . 

• a quantum effect which survives these averages is the behavior of the cross-section 
at very low energy. This is the quantum threshold regime, which extends up to an 
energy of the order of U /i~^/^Cg for an potential with a Cq/t^ long range, following 
Julienne and Mies J38| . This effect explains the very large helium-helium cross-section 
at low energy [10|. For argon-argon and lithium-argon collisions, we find that this 
regime could be observed for temperatures of the order of 10~^ and 2 x 10~^ Kelvin 
respectively, well below the lowest temperatures obtained in argon expansions, of the 
order of 1 Kelvin ll|, |l9| . 



B. Calculations with a Lennard- Jones 12 — 6 potential 

We consider an interaction potential of the 12 — 6 Lennard- Jones type: 

(37) 






n^^ i-Zi /v«0 



ai2 a^ 



ly 1-^ /v'O 



using the usual reduced quantities (see paragraph 8.2 of [ISi)- In a first step, we have 
calculated the deflection function x{b) as a function of the reduced energy E* = fig'^/(2e), 
by numerical integration. Our calculation reproduces well the analytic deflection function 
of BV for E* <0.1 (see appendix A of |ll|). 

We have then calculated the Q^^\g) cross-sections with / = 1 and / = 2 as a function of 
E* . In the low-energy range, i?* -C 1, these two cross-sections are well approximated by 
an [E*)~^/'^ behavior and, to make this behavior very clear, we have plotted in figure 5 the 
variation of Q^''\g)E*^'^ as a function of E*. At low energy, our calculation presents a small 
numerical noise associated to the orbiting singularity. We have not tried to reduce this noise 
by introducing some analytic approximations of the function x(6) near this singularity. 

When E* — > 0, we get: 

Q^'\g) ^ (7.81 ± 0.01)(t'E*-^/^ (38) 
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and 

g(2) {g) ^ (5.95 ± O.Ol)^^^*-!/^ ^39) 

Then, we have calculated numerically the fi*^''*-' integrals as a function of the reduced 
temperature T* = ksT/e and we have verified their low energy limits analytically (this is 
easy as the integral in equation (jHBI) can be expressed by a Gamma function). The various 
f2(''*)(T) take a form analogous to equation ((Tj): 

fi(2.2)(T) = c(''^)(T*) f^V ' ("t^V^' (40) 

\ /i y KkBTJ 

We have plotted in figure 6 the variations of three functions C^^'^'^T*). We are mostly 
interested in their low-temperature limits, which are necessary for our calculations. We 
thus obtain C^^'^^(O) = 3.00, in excellent agreement with the value 2.99 obtained by 
BV [11], C(i'2)(0) = 3.94 and C(i'^)(0) = 1.48. We also get the value of the ratio 
p„ = C(^'2)(0)/C(2-2)(o) = 1.312 ±0.004, also in very good agrement with the value po = 1.32 
deduced from the analytic deflection function of BV jllf . 



VIII. APPENDIX B: CALCULATION OF F 2 AND Fj 

We must calculate: 



ksUc 
with X = II and ±. We introduce the center of mass velocity Vc 



^x,2 ^ j^ J{g^^)dQ J AEx,2fi{^^i)f2{^2)d'^r,d\2dn (41) 



Vcm = avi + /5V2 (42) 

with a = rrii/M, (3 = m2/M and M = rrii + m2. The relative velocity before the collision 
is g = vi — V2 and after the collision is g/. The deflection angle noted x is the angle 
between g and g/, which have the same modulus. We must also introduce the azimuth cp 
measuring the orientation in space of the plane g, g/. This azimuth ip is equally probable 
in the < ip < 271 range. The final velocity of atom 2 is given by: 

V2/ = avi + I3v2 - agf (43) 

We express its parallel and perpendicular components as a function of the initial velocities, 
of the relative velocity and of the various angles and we want to calculate: 

AE2X = ^ (v2^,^ - V2\) (44) 

We then integrate these quantities over the angles and over the velocities vi and V2 with 
the normalized functions /i(vi)/2(v2). All the terms involve double products of the type 
VixVjYi where i,j stands for atom indices 1,2 and X^Y for || or _L. After integration, all 
these products vanish, excepted the ones with i = j and X = Y and their values are simply 
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related to the temperature associated to this atom and this degree of freedom: 



AR. = ^ 



kBTui 



ll^a^d 



nil 



cosxY 



kj^I ||2 / 2 2 n / 

—(a sin X + 2q;p (1 — cosx 

7712 ^ 



(kBT^l kBT^2\ 2 • 2 

+ I \ I a sm X 



nil 



ni2 



(45) 



AE2, = ^ 



m2 
2" 



A;bT||i , A;bT||2\ 2-2 , ^b^xi 2 /o , 2 . 

a sm xH a (-3 + cos x~4cosx 



+ 



nil 



ni2 



mi 



kuT 



BJ-±2 



ni2 



(Aaf3 (1 — cos x) + tt^ sin^ x 



(46) 



Finally we integrate over x- The integrals involve the angular weights (1 • 



and 2 only, i.e. they involve only Qi 2 and Qil defined by equation 



cos' x) with / = 1 
There was no 

terms involving / = 1 terms in the pure gas case. This new feature is due to the fact that 
the forward-backward symmetry existing in the pure gas case is now broken by the different 
masses of the two colliding atoms and by the fact that we calculate the energy transfer for 
atom 2 only. 

We have made an approximation by considering gQ 1^2(9) as independent of g and we must 
fix its value. To be coherent with the pure gas case, we have made the same approximation 
in the pure gas case and we have identified the two values of A{z): 



){i) 



32. 



)(',2), 
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(47) 



As fii 2 (Tm) varies very slowly with T^, like T^^, the exact choice for T^ is not important 
and we have always used Tm = (T\\i + 2Tj_i)/3. 
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FIG. 1: The reduced temperatures tx-i are plotted as a function of the reduced z-coordinate C- For 
Tim and r^i, this plot is identical to figure 5 of Beijerinck and Verster [5]. The ultra cooling appears 
on r||2, which reaches a terminal value lower than tim . The parallel and perpendicular temperatures 
of the seeded gas separate for larger Q values than the temperatures of the carrier gas and this fact 
is the basis of the ultra cooling of the seeded gas. The present calculation corresponds to the case 
of lithium seeded in argon, with ps = 2.55 and po = 1.32, and it predicts 2]|2oo/2]|ioo = 0.38. 

FIG. 2: The calculated parallel temperature ratio T\\2co/T\\i^ is plotted as a function of the mass 
ratio 1712 /ttii, spanning the range from hydrogen to rubidium seeded in argon. Each curve corre- 
sponds to a different value of the ratio ps and is labeled by this value. 



FIG. 4: The parallel temperature of lithium seeded in argon as a function of argon pressure in 
the source. The points represent our measurements. The full curve is the result of our theoretical 
calculation with Tu2/Tui = 0.38. The agreement between theory and experiment is quite good. The 
dashed curve results from a fit of the four data points with an argon pressure below 400 millibar: the 
best fit corresponds to a temperature ratio Tj|2/Tj|i = 0.31. For larger source pressures, the lithium 
parallel temperature remains roughly constant, in disagreement with theory (see the discussion in 
the text). 



FIG. 5: Plot of the variations of y^'^ = Q^^\g)E*^/'^ /a"^ as a function of ii^* for / = 1 (upper curve) 
and / = 2 (lower curve). Both quantities approach their limit rapidly when E* < 0.5. 



FIG. 6: Plot of the variations of the quantities (7' '^^(T*) defined by equation ()4U() as a function 
of T* = /c^T/e for (/,s) = (1,1) (dots); (1,2) (dash-dotted curve) and (2,2) (solid curve). 



FIG. 3: Laser induced fluorescence signal as a function of the laser frequency: the dots represent 
the experimental data and the full curves are the Gaussian flts to each line. The Fabry-Perot 
used for calibration has a free spectral range equal to 251.4 it 0.5 MHz. The saturated absorption 
spectrum providing lithium Doppler-free peaks is represented for frequencies above 1250 MHz. 
The fluorescence peaks are labeled by the ground state F value (F = 1 or 2) and by a letter 
corresponding to the laser beam. The angle between the atomic beam and the laser beam is 
61a = 47.9 ± 0.5° for beam A and Ob ~ 90° for beam B. 
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